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Abstract — Concurrent synchronization is a regime where di- 
verse groups of fully synchronized dynamic systems stably coex- 
ist. We study global exponential synchronization and concurrent 
synchronization in the context of Lagrangian systems control. 
In a network constructed by adding diffusive couplings to 
robot manipulators or mobile robots, a decentralized tracking 
control law globally exponentially synchronizes an arbitrary 
number of robots, and represents a generalization of the average 
consensus problem. Exact nonlinear stability guarantees and 
synchronization conditions are derived by contraction analysis. 
The proposed decentralized strategy is further extended to 
adaptive synchronization and partial-state coupling. 



I. Introduction 

Distributed and decentralized synchronization phenomena 
of large dynamic groups are areas of intense research. In this 
article, synchronization is defined as a complete match of all 
configuration variables of each dynamic system such that xi = 
X2 = • • • = x p and p denotes the number of sub-systems in the 
network. While such a definition directly concerns the angular 
position synchronization problem, this paper also addresses 
the synchronization of biased variables with application to the 
translational position coordination problem. In the latter case, 
synchronization corresponds to yi = y2 = • • • = y p where y$, 
1 < i < p connotes a vector of biased variables constructed 
from the configuration vector x$ such that y*(t) = x$(t) + 
hi(t) and the separation vector bi(t) is independent of the 
dynamics. Furthermore, we construct complex robot networks 
where multiple groups of fully synchronized elements coexist. 
Such concurrent synchronization seems pervasive in biology, 
and in particular in the brain where multiple rhythms coexist 
and neurons can exhibit many qualitatively different types of 
oscillations l32l . 

The objective of this paper is to establish a unified syn- 
chronization framework that can achieve both the synchro- 
nization of the configuration variables of the robots and the 
stable tracking of a common desired trajectory. Although 
an uncoupled trajectory tracking control law, in the absence 
of external disturbances, would achieve synchronization to a 
common desired trajectory, non-identical disturbances justify 
the mutual synchronization of the system variables. On the 
other hand, the synchronization to the average of initial 
conditions is not sufficient for multi-robot or multi-vehicle 
systems where a desired trajectory is explicitly given. For 
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example, a large swarm of robots can first synchronize their 
attitudes and positions to form a certain formation pattern, 
then track the common desired trajectory to accomplish the 
given mission. In production processes, such as manufacturing 
and automotive applications, where high flexibility, manipu- 
lability, and maneuverability cannot be achieved by a single 
system l36lL there has been widespread interest in cooperative 
schemes for multiple robot manipulators that track a prede- 
fined trajectory. A stellar formation flight interferometer (8), 
l9l is another example where precision control of relative 
spacecraft motions is indispensable. The proposed synchro- 
nization tracking control law can be implemented for such 
purposes, where a common desired trajectory can be explicitly 
given. The proposed strategy can achieve more efficient and 
robust performance through local interactions, especially in 
the presence of non-identical external disturbances. Further, 
we generalize the proposed control law such that multiple 
dynamic systems can synchronize themselves from arbitrary 
initial conditions without the need for a common reference 
trajectory. As a result, other potential applications include 
oscillation synchronization of robotic locomotion l33lL l37lL 
and tele-manipulation of robots 0, l27l . 

The main contributions of this work can be stated as follows. 

• Concurrent synchronization that exploits the multiple 
time scale behaviors from two types of inputs (a reference 
trajectory and local couplings) permits construction of a 
complex time-varying network comprised of numerous 
heterogeneous systems. 

• In contrast with prior work on consensus and flocking 
problems using graphs, the proposed strategy primarily 
deals with dynamic networks consisting of nonlinear 
time- varying dynamics. 

• We use contraction analysis l24l . l45l as our main nonlin- 
ear stability tool, thereby deriving exact and global results 
with exponential convergence, as opposed to asymptotic 
convergence of prior work. 

• The proposed control laws are of a decentralized form 
requiring only local velocity/position coupling feedback 
for global exponential convergence, thereby facilitating 
implementation in real systems. 

• The theory is generalized and extended to multi-robot 
systems with non-identical dynamics, linear coupling 
control, partial state coupling, uni-directional coupling, 
and adaptive control. 

A. Comparison with Related Work 

The consensus problems on graph l28l and the coordination 
of multi-agent systems (BJ, (23, EH, El are closely 
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related with the synchronization problem. In particular, the 
use of graph theory and Laplacian produced many interesting 
results C2I, |20|, (23, El, El, (H, 03. However, the 
synchronization to the average of initial conditions might 
not be directly applicable to multi-robot and multi-vehicle 
systems, where a desired trajectory is explicitly defined. A 
recent work l34l studied the consensus problem with a time- 
varying reference state, based on a single integrator model. In 
essence, the aforementioned work mainly deals with simple 
dynamic models such as linear systems and single or double 
integrator models without nonlinearly coupled inertia matrices. 
In contrast, we aim at addressing highly nonlinear systems 
(e.g. helicopters, attitude dynamics of spacecraft, walking 
robots, and manipulator robots). As shall be seen later, the 
proof of the synchronization for network systems that possess 
nonlinearly coupled inertia matrices is more involved. This 
paper focuses on such dynamic networks consisting of highly 
nonlinear systems. 

One notable work l36l introduced a nonlinear tracking 
control law that synchronized multiple robot manipulators in 
order to track a common desired trajectory. Due to its all-to-all 
coupling requirement, the number of variables to be estimated 
increases with the number of robots to be synchronized, which 
imposes a significant communication burden. Additionally, the 
feedback of estimated acceleration errors requires unnecessary 
information and complexity. Thus, a method to eliminate 
both the all-to-all coupling and the feedback of the accelera- 
tion terms is explored in this paper. Another nice approach 
to the synchronization of robot networks is to exploit the 
passivity of the input-output dynamics Q, (4). Its property 
of robustness to time delays is particularly attractive, while 
robot dynamics are passive only with velocity outputs unless 
composite variables are employed. In addition, the mutual 
synchronization problem, which not only synchronizes the 
sub-members but also enforces them to follow a common 
reference trajectory, is not addressed. In particular, it shall 
be shown that our proposed control law generalizes the robot 
control law presented in ID, while our convergence result is 
stronger (exponential). One recent work ifTHl used the passivity 
property for the path-following system to synchronize the path 
variables. The passive decomposition fT9lL l20l describes a 
strategy of decoupling into the internal group formation shape 
and the total group maneuver. Due to its dependency on a 
centralized control architecture, the decoupling is not generally 
ensured under the decentralized control. lf20l considers linear 
double integrator models. Another work ifTTTl proposed a nice 
control framework for controlling and coordinating a group 
of nonholonomic mobile robots, although the exact synchro- 
nization of individual vehicles was not pursued. One recent 
work 1 3 1 1 presented a framework called particles with coupled 
oscillator dynamics so that collective spatial patterns emerge. 
Our recent work also discussed more complex spacecraft 
dynamics to achieve the phase synchronization of spiral or 
circular translational trajectories in three dimensions as well 
as the synchronization of rotational dynamics. Also, it is 
important to note that the adaptive synchronization of multiple 
robotic manipulators in a single local coupling configuration 
was studied in (43J, l44l . 



The proposed synchronization framework using contraction 
analysis, first reported in (6), has a clear advantage in its broad 
applications to a larger class of identical or nonidentical non- 
linear systems even with complex coupling geometry including 
uni-directional couplings and partial degrees-of-freedom cou- 
plings, non-passive input-output, and time delays Q, |46l , 
while ensuring a simple decentralized coupling control law 
(see Figure [T] for network structures permitted here). The 
proofs and new results from [6] are expanded in this paper. 
In particular, concurrent synchronization that exploits two 
different types of inputs has not been studied in the literature. 

B. Organization 

Section [II] describes modeling of robots based on the 
Lagrangian formulation, and summarizes the key stability 



theorems. While Section III presents the main control law 
and its tracking stability, the proof of exponential synchro- 
nization is more involved and treated separately in Section 
IV The remainder of the paper further highlights the unique 
contributions of this work. Section IVl elucidates the concurrent 
synchronization of complex networks and the leader-follower 
problem. The main idea of this paper is extended to linear 
Proportional-Derivative (PD) couplings, and limited partial- 



state couplings in Section VI Key simulation results are also 



presented in Section VI 



II. Modeling and Nonlinear Stability Tools for 
Multi-Robot Networks 

A. Lagrangian Systems 

This paper is devoted to the use of the Lagrangian for- 
mulation for its simplicity in dealing with complex systems 
involving multiple dynamics. The equations of motion for 
a robot with multiple joints (q^ G W 1 ) can be derived by 
exploiting the Euler-Lagrange equations: 

Li = ~q/ Mi(qi)qi-Vi, — — - =7* 



dt dqi 



(1) 



where i, (1 < i < p) denotes the index of robots or dynamic 
systems comprising a network, and p is the total number of 
the individual elements. Equation ([I]) can be represented as 



Mi(qi)qi + Ci(qi, q^q* + g*(q*) n 



(2) 



where g;(q;) = ^r, and, ti is a generalized force or torque 
acting on the i-th robot. 

Note that we define C^q^q^) such that (M^ — 2C^) is 
skew- symmetric (38), and this property plays a central role in 
our stability analysis using contraction theory (7). 

The following key assumptions are used throughout this 
paper. Since the main applications of the present paper include 
fully actuated robot manipulators and spacecraft |9 ], the robot 
system in ([2]) is fully actuated. In other words, the number of 
control inputs is equal to the dimension of their configuration 
manifold (= n). The mass-inertia matrix M(q) is assumed to 
be uniformly positive definite, for all positions q in the robot 
workspace l38l . 
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B. Contraction Analysis for Global and Exponential Stability 

(1) Modular Stability Analysis: Although one popular method 
for modular stability analysis is to exploit the passivity formal- 
ism (38), we use contraction theory (24), (32), (39), (45) as 
an alternative tool for analyzing modular stability of coupled 
nonlinear systems. In particular, contraction analysis has more 
general and intuitive combination properties (e.g., hierarchies) 
than the passivity method, since it involves a state- space rather 
than an input-output method. 

(2) Differential State-State Analysis: Lyapunov's linearization 
method indicates that the local stability of the nonlinear system 
can be analyzed using its differential approximation. What is 
new in contraction theory is that a differential stability analysis 
can be made exact, thereby yielding global results on the time- 
varying nonlinear system. 

(3) Stronger Stability Results (Exponential and Global): For 
a robot dynamic model in ^ and a time-varying tracking 
control law, the straightforward use of Barbalat's lemma or 
LaSalle-Yoshizawa theorem (TS) yields asymptotic conver- 
gence results. While global exponential stability can be proven 
by using a Lyapunov function with a cross-term and additional 
constraints (T), Q, ifTTl . such a method is ad hoc, as compared 
to contraction theory. While exact exponential convergence 
might not be achievable in real systems due to the modeling 
errors, we believe that finding an explicit convergence rate 
with exponential stability, if possible, is important due to its 
superior tracking performance and property of robustness with 
respect to perturbations (e.g., see p. 339-350 in ifTTl ). 

A brief review of the results from (24) . (39), (45) is 
presented in this section. Note that contraction theory is a 
generalization of the classical Krasovskii's theorem (38) . and 
that approaches closely related to contraction, although not 
based on differential analysis, can be traced back to (Toll , (T3) 
and even to (22) . 

Consider a smooth nonlinear system 

x(t) =f(x(t),u(x,t),t) (3) 

where x(t) G R n , and f : R n x R m x R + — > R n . A virtual 
displacement, 5x is defined as an infinitesimal displacement 
at a fixed time- a common supposition in the calculus of 
variations. 

Theorem 1: For the system in if there exists a uniformly 
positive definite metric, 

M(x,t) = 0(x,t) T 0(x,t) (4) 

where is some smooth coordinate transformation of the 
virtual displacement, Sz = 0Jx, such that the associated 
generalized Jacobian, F is uniformly negative definite, i.e., 
3 A > such that 

F = (©(x,t) + ©(*,*) |Q ®(x,t) _1 < -AI, (5) 

then all system trajectories converge globally to a single tra- 
jectory exponentially fast regardless of the initial conditions, 
with a global exponential convergence rate of the largest 
eigenvalues of the symmetric part of F. 



Such a system is said to be contracting. The proof is given in 
(24). Equivalently, the system is contracting if 3 A > such 
that 

M+f — ) M + M— < -2AM (6) 
\ox J ax 

Equation ^ is useful for the stability proof of a Lagrangian 
system, since the inertia matrix M(q) of the robot dynamics 
in ^ can be chosen as the metric M in ([6|. 

The key advantage of contraction analysis is its superior 
combination property as follows. 

C. Contraction of Coupled Systems 

The following theorems are used to derive stability and 
synchronization of coupled Lagrangian systems. 

Theorem 2: Hierarchical combination 09l/ . Consider two 
contracting systems, of possibly different dimensions and met- 
rics, and connect them in series, leading to a smooth virtual 
dynamics of the form 

d_ fSzA = /F n \ (5zA 
dt \Sz 2 J \F 21 F 22 J \5z 2 ) 

Then, the combined system is contracting if F 2 i is bounded. 
Proof: see (7), ED, ED- ■ 

Theorem 3: Partial contraction i?5l/ . Consider a nonlinear 
system of the form x = f (x, x, t) and assume that the auxiliary 
system y = f(y,x,£) is contracting with respect to y. If a 
particular solution of the auxiliary y -system verifies a specific 
smooth property, then all trajectories of the original x-system 
verify this property exponentially. The original system is said 
to be partially contracting. 

Proof: See [45] for the virtual observer-like y system. ■ 

Theorem 4: Synchronization K45\l . Consider two coupled 
systems. If the dynamics equations verify 

xi - f(xi,t) = x 2 - f(x 2 ,t) 

where the function f (x, t) is contracting in an input- 
independent metric, then xi and x 2 will converge to each 
other exponentially, regardless of the initial conditions. Math- 
ematically, stable concurrent synchronization corresponds to 
convergence to a flow -invariant linear subspace of the global 
state space H32\l . 

Proof: This can be proven by constructing the virtual 
system y — f (y, t) = u(t) and using Theorem [3] (see (45) ). ■ 
Remark I: Whereas Theorem [2] can be proven for different 
metrics, Theorem [4] requires the same metric (e.g., inertia 
matrix) among the coupled systems. Hence, Theorem]?] cannot 
be directly applied to coupled Lagrangian systems. This is 
one of the motivations of the current paper, and elucidated in 
Section HV1 

III. Control Law and Its Tracking Stability 

A tracking controller introduced in this section achieves not 
only global exponential synchronization of the configuration 
variables, but also global exponential convergence to the 
desired trajectory. 
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(a) with bi-directional couplings 



A 



5-4 



^ ^ .IMA; 

(b) heterogeneous robots 



While the common desired time- varying trajectory (or the 
virtual leader dynamics) is denoted by qd(£), the reference 
velocity vector, q^ ?r is given by shifting the common desired 
velocity with the position error: 



= q d - Aqi = q d - A(q* - q d ) 



(9) 



(c) with uni- directional couplings (d) with uni-directional and bi- 
directional couplings 




A 



where A is a positive diagonal matrix. 

In contrast with l36l . the proposed control law requires only 
the local coupling feedback of the most adjacent robots (i — 1 
and i + 1) for exponential convergence (see Figure [I]). Note 
that the last (p-th) robot is connected with the first robot to 
form a ring network as suggested in l45l . Moreover, estimates 
of q are no longer required. 

The closed-loop dynamics using ^ and ([5]) become 

M(q i )s i + C(q i ,q i )s i + KiS i -K 2 s i _i-K 2 s i+ i =0 (10) 



(e) unbalanced graphs with feedback hierarchies (f) inline configuration where denotes the composite variable = q^ — q^ r 



Fig. 1. Network structures permitted in this paper. Networks in (a-d) are 
on balanced graphs, and each element has the same number of neighbors 
(i.e, regular graph). More complex structures, such as the unbalanced graphs 
shown in (e), can be constructed by concurrent synchronization presented in 
Section|V| The solid lines indicate the local couplings, whereas the dash lines 
indicate the reference input commands. 



A. Proposed Synchronization Control Strategy 

Let us first consider the robot networks shown in Fig. [TJa- 
d). The following control law, adopted from the single robot 
control law in l38lL is proposed for the i-th robot in a network 
comprised of p robots (p > 3). 

n =M(q i )q i>r + C(q i? q*)qi, r + g(q») (7) 

-Ki(t)(qi-q ijr )+ V — K 2 (t)(qj - q j?r ) 

which also permits the uni-directional couplings shown in 
Figures [T] (c) and (d). Note that m, which is the same for each 
robot, is the number of its neighbors that send the coupling 
signals to the i-th member, and the constant or time- varying 
set Mi(t) consists of such neighbors. This generalized form 
will be discussed again in Section [V] 

Equation ([7]) reduces to the following tracking control 
law for a two-way-ring symmetric structure, as shown in 
Figure [TJa) 

Ti =M(q i )q i , r + C(q i ,q i )q i , r + g(q*) - Ki(t)(q* - q*, r ) 
+ K 2 (t)(q i _i - qi_i, r ) + K 2 (t)(qi+i - qi+i,r) (8) 

where a uniformly positive definite matrix Ki(t) G R nxn 
is a feedback gain for the i-th robot, and another uniformly 
positive definite matrix K 2 (t) G R nxn is a coupling gain with 
the adjacent members (i — 1 and i + 1). While ([8} constructs 
a closed graph (e.g., ring), an inline configuration is also 
permitted since the first and last robot can simply connect the 
K 2 (£) gain from itself (see Fig. |TJf) and Section V-A). The 
above control law can also be applied to a network consisting 
of p non-identical robots (Figure [TJb)), as shall be seen in 
Section IV-Bl 



B. Modified Laplacian 

Let us define the following p x p block square matrices: 



j a,bJ 



A B 
B A 



BAB 
B A 



pxp 



A A 
A A 



A A 



pxp 



For a ring structure, defined from the controller in (|8j), [L^ B ] 
has only three nonzero matrix elements in each row (i.e., 
A,B,B). Then, we can write the closed-loop dynamics in 
([TO]) in the following block matrix form 

[M]x+ [C]x+ ([L^ i5 _ K2 ] + [U£ a ]) x = [U^Jx (11) 

where [M] = diag (M(qi), • • • , M(q p )), 

[C] = diag (C(qi, qi), • • • , C(q p , q p )), x = (sf , • • • , sj) . 

Definition [L^ _ K2 ] can be viewed as the modified Lapla- 
cian of the network in the context of graph theory. In other 
words, [L|^ _ K J indicates the connectivity with adjacent 
systems as well as the strength of the coupling by K 2 . Note 
that [Lj^ _ K J can be time- varying due to time- varying Ki 
and K 2 , or due to the switching topology, Ni(t) in 

Remark 2: The network graphs illustrated in Figure [TJa- 
d) are balanced since the in-degree of each node is equal 
to the out-degree f28l. The additional requirement for the 
stability analysis in this section is that the robots should 
be on a regular graph, where each member has the same 
number of neighbors. While this paper permits popular local 
coupling configurations (9), 031 from regular graphs, we will 
show that the assumption of a regular balanced graph can be 
relaxed in Section |V| In particular, unbalanced or non-regular 
graphs due to feedback hierarchies, as shown in Figure [TJ e )> 
can be employed by concurrent synchronization discussed in 
Section [VJ 

Remark 3: It should be noted that the matrix [L^ _ K J 
is different from the standard Laplacian found in |28| . By 
definition, every row sum of the Laplacian matrix on a 
balanced graph is zero. Hence, such Laplacian matrix always 
has a zero eigenvalue corresponding to a right eigenvector, 
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(1,1,- •• , 1) T . In contrast, a strictly positive definite 
is required for exponential tracking convergence 



-k 2 



for the proposed control law in this paper. In other words, 
unless otherwise noted, [L^ _ K ] is assumed to have no zero 
eigenvalue. For example, the block matrix for p = 4 becomes 



j k 1? -k 2 



+Ki -K 2 -K 2 

-K 2 +Ki -K 2 

-K 2 +Ki -K 2 

-K 2 -K 2 +Ki 



(12) 



which is positive definite for Ki 
is also true Vp, p > 3. 



2K2 > 0. This condition 



IV. Synchronization With/Without Tracking 

We prove the exponential synchronization of multiple La- 
grangian systems in this section. First, we describe the diffi- 
culties inherent in proving the synchronization of Lagrangian 



systems in Section |IV-A We then present the synchronization 



proof in Sections IV-C We also show that our method is more 
general than prior work by reducing our control law to the 
standard synchronization problem without trajectory tracking 



in Section IV-D| The adaptive synchronization is presented in 
Section H^E] 



C. Tracking Stability Analysis 

The following condition should be true for exponential 
tracking convergence to the common desired trajectory q^(t). 

Theorem 5: If [L^ K J is uniformly positive definite: 

[L^ 1) _ K J>0, Vt (13) 

then every robot follows the desired trajectory q_d(t) exponen- 
tially fast from any initial condition. 



Proof: We can cancel out the [U^ ] matrix term in (11) 
to obtain 



[M]x+[C]x+[L^_ K2 ]x = 0. 



(14) 



Equation ( [T4| ) corresponds to a conventional tracking problem 
with a block diagonal matrix of nonlinearly coupled inertia 
matrices, [M]. We use contraction theory to prove that x 
tends to zero exponentially and globally with [L^ _ K ] > 0. 
Consider the virtual system of y obtained by replacing x with 
y in ([14). 



[M]y+[C]y+[L^_ K2 ]y = 



(15) 



This virtual y system has two particular solutions: x 



and 0. The squared-length analysis with respect 



to the positive-definite metric [M] yields 

d 



dt 



(<5y T [M],5y) = 28y T [M]8y + Sy T [M]Sy (16) 



= -2Sy T {[C]Sy + \L p Ku _ K2 ]6y) + 5y T [M]Sy 
= -28y T [L^_ K2 ]5y 

where we used the skew-symmetric property of [M] — 2 [C] . 

Accordingly, [L^ _ K ] > will make the system contract- 
ing (5y — > 0), thus all solutions of y converge to a single 
trajectory globally and exponentially fast (Theorems [T] and 
[3J. This in turn indicates that the composite variable of each 
robot tends to zero exponentially (s$ — > 0). By definition of 
Si = — qd + A(qi — q^), this show global exponential 
convergence of q^ to the common reference trajectory q^(t) 
(see also Theorem [2]). ■ 

Given Ki > 0, K 2 > 0, it can be shown that a sufficient 
condition for the positive-definiteness of [L^ _ K ] is Ki — 
K 2 > for p = 2, and Ki - 2K 2 for p > 3. 1 ' 

The next question to be addressed is how to guarantee the 
synchronization of the individual dynamics. 



A. Challenges with Nonlinear Inertia Matrix 

The difficulties associated with nonlinear time-varying in- 
ertia matrices can be easily demonstrated with the following 
two-robot example. The closed-loop dynamics of two identical 
robots from ( [TO] ) becomes 

M(qi)si + C(qi, qi)si + (K x + K 2 )s x = u(t) 
M(q 2 )s 2 + C(q 2 , q 2 )s 2 + (Ki + K 2 )s 2 = u(t) 



(17) 



Note that u(t) = K 2 (si + s 2 ) and that s$, i = 1,2 is the 
composite variable defined in ( [TQ| ). 

Direct application of synchronization (Theorem]?]) appears 
elusive since we have to prove that (17) are contracting in 
the same metric while preserving the input symmetry l32ll . 
Hence, this can be viewed as a higher order contraction 
problem (25). For example, multiplying (17) by M _1 breaks 
the input symmetry: i.e., M -1 (qi)u((:) ^ M _1 (q 2 )u( ( :). In 
essence, M(qi) / M(q 2 ) makes this problem intractable in 
general. 

Instead, assume that M(q) becomes a constant matrix, 
thereby making C(q, q) zero. Then, we can easily prove Si 
and s 2 tend to each other from 

Msi + (Ki + K 2 )si = u(t), Ms 2 + (Ki + K 2 )s 2 = u(t) 

(18) 

The following virtual y-system with the common input u(t) 



My + (K 1 +K 2 )y = u(t) 



(19) 



is partially contracting with Ki + K 2 > (see Theorem [3J. 
Hence, its particular solutions si and s 2 tend to each other 
exponentially fast according to the synchronization theorem 
(Theorem]?]). Without loss of generality, this result can easily 
be extended to arbitrarily large networks. The synchronization 
of a large network with a constant metric, as seen in ( [T8] ), is 
already discussed in l45l using contraction analysis. 

We now turn to a more difficult problem focused on the 
synchronization of two robots with non-constant nonlinear 
metrics (M(qi) ^ M(q 2 )). 

B. Contraction with Multiple Time Scales 

In this section, we show that we can render the system 
synchronized first, then follow the common trajectory by 
tuning the gains properly. This indicates that there exist two 
different time scales in the closed-loop systems constructed 
with the proposed controllers. This multi-time-scale behavior 
will be exploited in the subsequent sections. 
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Recall the closed-loop dynamics given in ( [T4| ). Since 
_ K ] is a real symmetric matrix, we can perform the 
spectral decomposition [5]. This is a special case of the con- 
current synchronization [32] that corresponds to convergence 
to a flow invariant subspace (the eigenspace). 

[L^,-kJ=V[D]V t (20) 

where [D] is a block diagonal matrix, and the square matrix 
V is composed of the orthonormal eigenvectors such that 
V T V = VV T = I pn , since the symmetry of [L^ i _ K J gives 
rise to real eigenvalues and orthonormal eigenvectors [42]. 
Pre-multiplying ( [l4| ) by V T and setting V T x = z result in 

(V T [M]V) z + (V T [C]V) z + [D]z = (21) 

Then, we can develop the squared-length analysis, as in 



Sj=S 2 synchronization 



tracking 




(16). This follows from the fact that (V T [M]V) is always 



symmetric positive definite due to a symmetric positive defi- 
nite [M]. 

Since the modified Laplacian [L^ _ K ] represents a reg- 
ular graph, where each member has the same number of 
neighbors (= 2 for p > 3), 

[1] = -^L[I n ,I n ,... ,I n ] T (22) 

VP 

is the pnx n block column matrix of eigenvectors associated 
with the eigenvalues A(Ki — 2K2) for p > 3. Note that I n 
denotes the n x n identity matrix, and the [1] matrix consists 
of p matrices of I n . The eigenvector matrix [1] represents the 
common reference trajectory tracking state. 



Fig. 2. Multiple timescales of synchronization (faster) and tracking (slower). 
The dashed line indicates the desired trajectory, and arrows indicate increasing 
time. The drawing is conceptual, since strictly speaking si and S2 synchronize 
exponentially. 



C. Stability Analysis of Exponential Synchronization 

Using the results from the previous sections, we present the 
main theorem on synchronization. 

Theorem 6: Assume that the conditions in Theorem |3] are 
true, so that the individual dynamics are exponentially tracking 
the common desired trajectory. A swarm of p robots synchro- 
nizes exponentially from any initial conditions if 3 diagonal 
matrices Ki > 0, K 2 > 0, A>0 such that 

D 2 = V^ nc [L^_ K2 ]V sync >0, Vi (25) 

Proof: Consider the virtual system of y from pi) : 

(V T [M] V) y + (V T [C] V) y + [D]y = 0, (26) 

This has y = V T x and y = as particular solutions, which 



We can define a pn x (p - l)n matrix V sync that consists can be written in terms of ? = fat, vl 



of the orthonormal eigenvectors other than [1] such that 

[i] T 



V J V 



([1 

i] T [i] [i] T v, 



sync 



(23) 



yt = [i] T x 



V 7 x 

v sync 



and 



yt 
y s 



(27) 



V 1 [11 

. v sync l- 1 -} 



v 



sync 
V 

sync v sync 







(p-l)nxn 



nx (p—l)n 
t(p-l)n 



where we used the orthogonality between [1] and V sync . 

Hence, the block diagonal matrix [D] , which represents the 
eigenvalues of [L^ Ci _ K J, can be partitioned from (20) 

fDl =V T \U 



— K 

k„-kJV 

kJ[1 

-K 



. v sync l Ki 



][i] 



[l] T [L P Kl) - 

v sync l Ki , 



K 2 J 
K 2 



v sync 

]v. 



sync 







Or 



nx (p— l)n 

D 2 



(24) 



^(p-l)nXn 

It should be emphasized that Di, which equals Ki — 2K2 
for p > 3, represents the tracking gain, while D 2 corresponds 
to the synchronization gain. We can choose the diagonal 
control gain matrices Ki and K 2 such that 

D 2 = VlM u - K2 ]Vsync > Di = [l] T [L^ ii _ K J[l], 

thereby ensuring that the robots synchronize faster than they 
follow the common desired trajectory. 

This multi timescale behavior is graphically illustrated in 
Figure [2] The figure depicts that Si and s 2 synchronize first, 
then they converge to the desired trajectory while staying 
together. This observation motivates separation of the two 
different time scales, namely Di for tracking and D 2 for 
synchronization. 



For [D] > 0, we can show that the above virtual system is 
contracting (i.e., Sy — > globally and exponentially). We take 
the symmetric positive definite block matrix V T [M]V as our 
contraction metric. 

Performing the squared-length analysis with respect to this 
metric yields 

j t 8y T (V T [M]V) Sy = -2Sy T ((V T [C]V) Sy + \D]Sy 
+ Sy T (V T [M]V) Sy = -2Sy T [F>}Sy (28) 

where we used the skew-symmetric property of ^V T [M] V^j — 
2(V T [C]V). 

The above equation can be rewritten in terms of two 
different time scales 



d 
dt 



8y t 
5y s 



[1} T [M}[1] [1] T [M]V 



£y s 



sync 



Di 




[M][l] 


D 2 



Vf ync [M]V. 
Sy s 



sync 
sync 



$y s 



(29) 



If Di > and D 2 > 0, the combined virtual system in 
( [26] ) is contracting. In other words, Sy — > exponentially fast. 
This in turn implies that all solutions of y tend to the single 
trajectory. In particular, the tracking (Sy t — > 0) is associated 
with Di, and synchronization (Sy s — > 0) is associated with 
D 2 . 
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s p ) and VjL^x from 



sync 

s„ — > has 



As a result, [l] T x = l/y/p(si + • 
( [27] ) tend to zero exponentially. Note that si 
already been proven with Di > for Theorem |5J which is a 
sufficient condition to make the sum of the composite variables 
also tend to zero (i.e., [l] T x — > 0). What is new in this section 
is the synchronization V^ nc x — > and its convergence rate 
with D 2 > 0. 

It is straightforward to show that Vf ync x — > and A > 
also hierarchically make qi , • • •> Qp synchronize globally 
exponentially fast (see Theorem [2]). This can be verified by 
the following contracting dynamics constructed from ( [T0| ) 



sync ) 



sync 



{q} 



V 1 x 

sync 



> 
(30) 

where {q} = (qf , • • • , qJ) T . Note that the orthonormal 
vectors V ' sync canceled the common input term q^ + Aq^. 
Also, [A] is a block diagonal matrix of A > 0, thereby yielding 
[A][l] =0 from 



v sync 



V syncV sync {c\} 



[l][l] T {q} = {q} 



(31) 



implies the synchronization 

1 = 1, • • • ,p. This also implies 



Consequently, Vj ync {q} - 
of the original state variable q^ 
that the diagonal terms of the metric, V^ nc [M] [1], tend to 
zero exponentially, thereby eliminating the coupling of the 
inertia term V T [M]V in (28 -29). This completes the proof 
of Theorem [6] ■ 

Remark 4: We assume here that qd(t) is identical for each 
member. If q_d(t) were different for each dynamics, — » Sj 
would imply the synchronization of q^ — q^ to the difference 
of the desired trajectories, which would be useful to construct 
phase synchronization of oscillatory trajectories (e.g., see 0, 
El). 

We can use Theorem [6] without finding V ' sync as follows. 



Corollary 1: The following condition, in lieu of ft25y, veri- 
fies Theorem^ 



K : 



J + [U^ 2 ]>0, V* 



(32) 



Proof: The block matrix [U^J also has [1] as its eigen- 
vector. We multiply [K sync ] = [L^ _ K J + [U^J by its 
orthonormal eigenvectors other than [1]: 



[K S2/nc ] sync — [-^K^ — K 2 ]^" s 



^ sync^-^2 



(33) 



which shows that V sync also represents the orthonormal 
eigenvectors of [K sync ] = [Lj^ _ K J + [U^-J. In other 
words, D2 corresponds to the eigenvalues of [K S2/nc ]. The 
remaining eigenvalue of [K sync ] that is associated with [1] is 
Ki + (p — 2)K2, which is greater than the tracking eigenvalue 
Di = Ki — 2K2 for p > 3. Hence, the synchronization occurs 
[L^.kJ + [U^J > 0. ■ 
Theorem [6] and Corollary [T] correspond to synchronization 
with stable tracking. As shall be seen in Section |V-B| multiple 
dynamics need not be identical to achieve stable synchroniza- 
tion. It is useful to note that the above condition corresponds 
to Ki + K 2 > for two-robot and three-robot networks 
(P = 2,3). 

Remark 5: Robust Synchronization. By extending l24lk we 
can show that coupled contracting dynamics have the prop- 
erty of robustness to bounded disturbances (7). It should be 



emphasized that exponential stability of contraction analysis 
facilitates such a perturbation analysis. In general, the proof 
of robustness with asymptotic convergence is more involved, 
or even leads to instability with respect to a certain class of 
perturbations (see p. 350 in 11711 ). 

In the next section, we show that a network of multiple 
robots can synchronize even without stable tracking. 

D. Synchronization Without a Common Reference Trajectory 

We now turn into the more standard synchronization prob- 
lem where the tracking gain is zero (Di = Ki — 2K2 = for 
P > 3), which fails the exponential tracking stability condition 
in Theorem |5j In this case, the modified Laplacian [L^ _ K ] 
reduces to the standard weighted Laplacian whose row sum is 
zero. For synchronization to the weighted average of the initial 
conditions, we do not require the common desired trajectory 
qd, and q^ can simply be set to zero as follows: 



q i>r = -Aqi 



Si = qi + Aq^ 



(34) 



In other words, our control strategy represents a more 
generalized framework for the synchronization of multiple 
Lagrangian systems. 

Theorem 7: Suppose that the individual tracking dynamics 
are indifferent. Hence, the conditions in Theorem [3] are not 
true. Nevertheless, a swarm of p identical robots asymptot- 
ically synchronize from any initial conditions if 3 diagonal 
matrices Ki > 0, K2 such that 



-K; 



; ] + [U£J>0, vt 



(35) 



For indifferent tracking, a common desired trajectory q_d(t) is 
no longer required. 

Proof: The combined virtual system per se is then semi- 



contracting (see l24l ) since the squared-length analysis in (28 
|29| ) yields the negative semi-definite matrix: 

T 



V = -2 



5y s 





D 2 



5y s 



< 



(36) 



While 5y t , representing the tracking dynamics, remains in a 
finite ball due to Di = 0, Sy s tends to zero asymptotically 
due to D2 > 0. This result can be proven as follows. Note 
that V = —^SyJ~DoSy s . V is uniformly continuous since a 
bounded Sy s from (26) leads to a bounded V. Due to V < 0, 
the use of Barbalat's lemma [38 ] verifies that V — > as t — > 
00. This implies that Sy s tends to zero asymptotically fast. 
As a result, V^ nc x tends to zero asymptotically. From the 
hierarchical combination discussed in ( [30| ), this also implies 

v sync IMJ 

indifferent tracking 

This will eventually decouple the metric matrix with A > 0, 
as seen in (|29|), since V^ irjr [M][l] tends to zero simultane- 

M(cfc) 
• turns 



S y ncl ^j 0. This completes the proof of Theorem [7] with 



sync 



ously as Vj^7 c {q} — > 0. As a result, when M(q^ 
is sufficiently close to zero, the convergence of Sy s 
exponential. 

Remark 6: Fast Inhibition. The dynamics of a large net- 
work with semi-contracting stability as in ( |36| ) can be instan- 
taneously transformed to contracting dynamics by the addition 
of a single inhibitory coupling link. In other words, a single 
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inhibitory link will also make Sy t — > 0. For instance, we can 
add a single inhibitory link between two arbitrary elements a 
and b while we keep the rest of the elements the same: B31 



r a =M(q a )q ar + C(q a , q a )<W + g(q a ) 

- 2K 2 s a + K 2 s a _i + K 2 s a+ i - K(s a + s b ) 
r b =M(q 6 )q 6r + C(q 6 , q 6 )q6r + g(q&) 

- 2K 2 s 6 + K 2 s 6 _i + K 2 s 5+ i - K(s a + s b ) 



(37) 



where 2K 2 is substituted for Ki, and K > 0. 

Hence, we can straightforwardly show that [L^ K J is 
now strictly positive definite, in contrast with the original 
semi-contracting system. As a result, the closed-loop system 
is contracting, resulting in — > and q^ — > from ([34]). 



This fast inhibition is useful to rapidly destroy unwanted 
synchronized oscillation of robots. 



E. Adaptive Synchronization 

We present the adaptive version of the proposed control 
law that adapts to the unknown parametric uncertainties of 
the robot dynamic models. Consider the following adaptive 
control law fT6lL l38lL l39l that has the same local coupling 
structure as the proposed control law in ([5]): 

Ti = Mi(qi)qi,r + Ci(qi, q;)q;, r + gi(q») - KiS* (38) 
+ K 2 s^_i + K 2 s i+ i = Yiki - KiSi + K 2 Si_i + K 2 s i+i 

where s$ = q^ — q^ r . Also, the parameter estimate a^ for the 
i-th member is updated by the correlation integral 



-tyTs 



(39) 



where T is a symmetric positive definite matrix. Hence, the 
closed-loop system for a network comprised of p non-identical 
robots can be written as 



[M] 






[r- 1 ] 



X 

{a} 



[C] 




x 

{&} 



[ L K!,-K 2 



[Y]' 



{a} 



(40) 



= 



where [M] and [C] are the block diagonal matrices of 
Mi(qi) and C^q^q^), i = 1, • • • as defined in (11 



The additional block diagonal matrices are defined from 
(39} such that [I 1 " 1 ] = diag(r _1 ,r _1 • • • .T' 1 )^ \Y] = 



V T [L^ Ci _ K JV = [D], to d40j> leads to the following virtual 
system of (y?,y%) T EE M 



V T [M]V 





(41) 



= 



The virtual system has two particular solutions: 
: V T x, y 2 = {a}, and yi = 0, y 2 = 
virtual length analysis indicates that (PT 



yi = 

The virtual length analysis indicates mat (1411) is semi- 
contracting by the negative semi-definite Jacobian with [D] > 
0: 



T r 







Syi 

^2 



(42) 



dV = _ 2 (S yi 
dt \Sy 2 , 

Using Barbalat's lemma (see Section IV-D| ), it is straightfor- 
ward to show that Syi tends asymptotically to zero from any 
initial condition. Also, [D] can be decomposed to the tracking 
and synchronization gains, and the rest of the proof follows 
Theorems 6] and [7] Consequently, the adaptive synchronization 
law in ([38 ) synchronizes the states of multiple dynamics in the 
presence of parametric model uncertainties. ■ 
While the synchronization of the estimates of the physical 
parameters (Jy 2 — > 0) is not automatically guaranteed due to 
the semi-contracting stability, the additional condition of the 
persistency of excitation [ 38 ] leads to the convergence of a to 



zero. Results of the simulation are presented in Section VII-C 



F. Examples 

For the case of a two-robot network, we can easily verify 



Z = V T X : 

Dx = Ki 



X : 



41 

V2 l 



sync J 

K 2 , D 2 = Ki + K 2 



A 1 , 



(43) 



If Ki 
V in (IS 

V 



- K 2 > Ki - K 2 > 0, the rate of the virtual length 
is uniformly negative definite: 

T 



5y s 



Kx-K 2 






Kx+K 2 



fyt 
$y s 



< o. 



Consequently, the combined virtual system in ( [26] ) is contract- 
ing. As a result, si + s 2 — > and si — > s 2 exponentially. It 
is straightforward to show that si — > s 2 also hierarchically 



makes qi tend to q 2 exponentially as in (30). 

The results hold for arbitrarily large networks as well. For 
example, a network of three robots has the following V, whose 



diag(Yi, Y 2 , • • • , Y p ). Also, x = (sf , s% , ■ • • , sJ) T , and columns are orthonormal eigenvectors of [L; 
{a} = (af , a|\ • • • ,aJ) T where a^ denotes an error of the 



p=3 

K 1? -K 2 



parameter estimate such that a$ = a^ — a^. Note that a^ is a 
constant vector of the true parameter values for the i-th robot, 
resulting in a^ = a^. If each robot is identical, a^ = a for 

1 < i < p. 

Theorem 8: The adaptive synchronization law in pff] ) glob- 
ally asymptotically synchronizes the states of multiple dynam- 
ics in the presence of parametric model uncertainties if the 
condition dJ5] ) holds. 

Proof: Similar to Section IV-C| applying the spectral 



V3 A 



Vs V6 
V3 1 V6 1 



V2 1 - 

7I 1 J 



The block diagonal matrix [D] is also computed as 
diag (Ki - 2K 2 , Kx + K 2 , K x + K 2 ). 

For a four-robot network p = 4 (see Figure [IJa)), [D] = 
diag(Kx - 2K 2 ,Ki +2K 2 ,Ki,Ki), and 



transformation, using the augmented V a = diag (V, l pn ) and 



II -a/2 I 

I -I -y/2I o 

II y/21 
I -I y/21 



By inspecting the associated eigenvectors, we can notice that 
Ki represents the synchronization gain associated with the 
synchronization of diagonal members (qi = q3,q2 = q4), 
while Ki + 2K 2 represents the synchronization gain of direct 
couplings (e.g., qi = q 2 ,q3 = q4). This is a percolation 
effect discussed in (32). The percolation effect can be exploited 
in order to prove the synchronization to the large unknown 
invariant set (e.g., V^ nc x = from a large network) by 
verifying the synchronization to a known, not necessarily 
invariant, subset of the global flow-invariant set (T2l . 

V. Concurrent Synchronization of 
Heterogeneous Groups on Unbalanced Graphs 

We further generalize the proposed synchronization frame- 
work in the context of the concurrent synchronization of mul- 
tiple heterogeneous networks and leader-follower networks, 
which permit construction of complex network structures. 

A. Inline Configuration and Uni-Directional Digraph 

We can also consider the inline configuration shown in 
Fig. [TJf), when maintaining a ring structure is not feasible 
possibly due to communication problems. 

Corollary 2: All the previous theorems hold for network 
structures on inline configuration, if the first and last robot 
adjust the proposed control law in ([3]) as 

n = M(q i )q i5r + C(q i ,q i )q i , r + g(q i )-(Ki-K 2 )s i + K 2 s i 

(44) 

where i = 1 or i = p, and j is the index for its sole neighbor. 
Proof: We can show that the modified Laplacian 




2 nd Network 




/> 

* * x - 



Fig. 3. Concurrent synchronization between two different groups. The desired 
trajectory inputs are denoted by the dashed-lines whereas the solid lines 
indicate mutual diffusive couplings. The independent leader sends the same 
desired trajectory input to the first network group. If we view the dashed-lines 
as edges of the graph, the network is on an unbalanced graph. 



where M^, C^, and gi(q*) represent the i-th robot dynamics, 
which can be different from robot to robot. Each robot has the 
same number of configuration variables (q^ G W 1 ). 

Corollary 4: The proposed tracking and synchronization 
control law in ([3]) can easily be applied to a network consisting 
of heterogeneous robots in ([2]) if the stable tracking condition 
in Theorem\6\is true. 

Proof: The M and C matrix notations used in the 
previous sections can be interpreted as M(qi) — > Mi(qi) 
and M(q 2 ) M 2 (q 2 ) with Mi(-) ^ M 2 (-) (the same 
for the C matrices). Hence, the assumption of non-identical 
dynamics does not alter the proof of exponential tracking in 
Theorem [5] and stable synchronization in Theorem [6] However, 
the synchronization with indifferent tracking (Theorem [7]) is 



I L Ki,-k 2 ] constructed from both (|8]) and (|44| is symmetric no longer true for non-identical robots since Vf ync x = is 

not a flow-invariant manifold, and qi = q 2 does not cancel 
the off-diagonal terms 1/2 [Mi(qi) — M 2 (q 2 )] in the metric 
matrix V T [M]V. ■ 



and still has [1] as its eigenvector associatea its eigenvalue 
Di = Ki — 2K 2 . Then the proofs of the theorems easily 
follow. ■ 

Let us now consider the generalized control law in ([7]), 
which permits regular digraphs, as shown in Figures [TJc) and 
(d). The modified Laplacian matrix [L^ _ K ] defined from 
([7]) might have two or more than three nonzero elements in 
each row. Further, it is no longer symmetric. 

Corollary 3: All the previous theorems are valid for regular 
graphs with uni-directional couplings or both uni-directional 
and bi-directional couplings. 

Proof: It is straightforward to verify that the sym- 



metric matrix 



K 2 



[LL,-kJ T has the 



same 



number of nonzero elements at each row (i.e., reg- 
ular). Hence, we can conclude that all the previous 
proofs still hold. For example, the synchronization occurs 
when V^ nc ([L^ i5 _ K2 ] + [L^ i5 _ K2 ] r )v s , nc > while 

I 1 ] T ([ L k 1? -kJ + [ L k 15 -kJ T ) [1] determines the stability 
and convergence rate of the trajectory tracking. ■ 

B. Synchronization of Heterogeneous Robots 

Consider the proposed control law for a network comprised 
of heterogeneous dynamics as follows: 

Ti = M^q^q^ + Ci(qi, qi)q», r + g»(q») (45) 
- Ki(qi - q i>r ) + K 2 (qi_i - q;-i, r ) + K 2 (q i+i - q i+ i, r ) 



C. Concurrent Synchronization of Heterogeneous Groups 

We present a new method of achieving the concurrent 
synchronization of multiple heterogeneous networks in this 
section. Such a method can be used to apply the main 
results in the previous sections to more complex networks, as 
shown in Figure [TJe). The networks in the figure are neither 
regular nor balanced due to the reference input couplings. 
In l32lL concurrent synchronization is defined as a regime 
where the ensemble of dynamical elements is divided into 
multiple groups of fully synchronized elements, but elements 
from different groups are not necessarily synchronized and can 
exhibit entirely different dynamics. 

As discussed in the previous sections, we pay particular 
attention to the fact that there exist two different time scales 
of the proposed synchronization tracking control law. In 
particular, one rate is associated with the trajectory tracking 
(Di), while the other represents the convergence rate of the 
synchronization (D 2 ). This implies that there are two different 
inputs to the system, namely, the common reference trajectory 
q_d(t) and the diffusive couplings with the adjacent members. 
Accordingly, we exploit a desired trajectory qd(t) to create 
multiple combinations of different dynamics groups. 
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For instance, Figure [3] represents the concurrent synchro- 
nization of two different dynamic networks. The first network, 
consisting of four heterogeneous robots, has the diffusive 
coupling structure proposed by the tracking control law in ([SJ. 
The independent leader sends a desired trajectory command 
qd(t) to each member of the first network. Note that the 
dimension of the leader dynamics need not be equal to that 
of the first network, since a signal from the leader can be 
pre-conditioned. With an appropriate selection of gains, each 
dynamics in the first network synchronize while exponentially 
following the leader. Therefore, the proposed scheme can be 
interpreted in the context of the leader- follower problem (O, 

ED, El). 

The second network consists of three heterogeneous dynam- 
ics, also different from those of the first group. Each element 
receives a different desired trajectory input from the adjacent 
element of the first network. Again, the desired trajectory input 
for the second group can be pre-conditioned from the first 
network, resulting in 



q2d(*) = f(qi^) 



(46) 



where qi is the state vector of the member in the first network 
that sends in the desired command and f is a differentiable 
nonlinear function of qi . As a result, the dynamics of the two 
networks can be entirely different (i.e., qi G M m , q2 G M n , 
and m ^ n). Note that q 2 ^, needed for the proposed control 
law in ([7]), can be obtained from a nonlinear observer l38l . 
The concurrent synchronization between the first and second 
networks can be proven as follows. 

Corollary 5: Consider the network structures that are un- 
balanced (i.e., the in-degrees are not equal to the out-degrees) 
due to the directional reference input connections among 
the regular balanced graphs, shown as the dashed lines in 
Figures ^e) and [5] The robots on such graphs globally 
exponentially synchronize if the individual groups synchronize 
from Theorems ^5\and^ 

Proof: For the example shown in Figure [3j once the 
first network synchronizes, the second network also ends up 
receiving the same desired trajectory to follow, while they 
interact to synchronize exponentially fast. Accordingly, we 
can achieve concurrent synchronization between two different 
network groups. The proof of Theorem [6] holds until (30), 
where now the desired trajectory inputs are different for each 
robot. Then, we can conclude that the robots in the second 
network synchronize, if the desired reference inputs q 2 ^,i, 
1 < i < p, sent from the first network, synchronize: 



Vf y „ c {q 2 } + (Vf y „ c [A]V 

— (V syncl^l i ' ' ' i r p ] * 



sync. J V sync {q.2} 

= V T Xo 

v sync^* 



(47) 







where = q 2 d,i + Aq 2 d,i, and x 2 is the vector of the 
composite variables (s$) of the second network. If the first 
network synchronizes, the reference trajectory term vanishes 
(V 



T r T . 
sync l 1 i 



tion of the second network (i.e., V^ nc {q 2 } 



0), thereby ensuring the synchroniza- 
0). This can be 

extended to arbitrarily large groups of synchronized dynamics 
by appropriately assigning the desired trajectory inputs and 
the diffusive couplings. ■ 



Remark 7: The previous example in Figure [3] and Corol- 
lary [5] can be interpreted in terms of feedback hierarchies 
(see Theorem [2]). The first network in the figure provides 
feedforward commands to the second network, but does not 
receive any command from the second network. Note that 
the dynamics at each level can be very different, in order to 
construct a dynamic concurrent combination of heterogeneous 
networks. Since the feedforward inputs can be appropriately 
scaled and conditioned, the dimensions between hierarchical 
layers can also be very different. For example, the dynamics 
higher in the hierarchy need not be oscillators, but could 
be systems with multiple equilibria (e.g., x = — VV, with 
local minima in V (39)), with synchronization correspond- 
ing to convergence toward a common minimum. Concurrent 
synchronization discussed in this section can be exploited to 
construct a large complex network consisting of heterogeneous 
dynamics such as robots, ground vehicles, and unmanned 
aerial vehicles. 

VI. Extensions 
Let us further extend the proposed control law. 

A. Synchronization of Robots Using Linear PD Control 

One may consider the following Proportional and Derivative 
(PD) coupling control law for two identical robots from ^ 
with p = 2, 



Ti = -Ki(qi + Aqi) + K 2 (q 2 + Aq 2 ) 
r 2 = -Ki(q 2 + Aq 2 ) + K 2 (qi + Aqi) 



(48) 



where q* = q* — q^ and the bounded reference position q^ 
has zero velocity such that q^ = q^. 

For simplicity, the gravity term in ^ is assumed to be zero or 
canceled by a feed-forward control law. Then, the closed-loop 
dynamics satisfy 



M(qi)qi + C(qi, qi)qi + Kqi + KAqi = u(t) 
M(q 2 )q 2 + C(q 2 , q 2 )q 2 + Kq 2 + KAq 2 = u(t) 

where K = K x + K 2 and u(t) = K 2 (qi + q 2 ) + K 2 A(qi 



(49) 



Similar to Section IV-C we can perform a spectral decom- 
position: 

(V T [M] V) V T x + (V T [C] V) V T x (50) 
+ (V T [L^ K15 _ K2 ]V) V-x+ (V-[L^ K1A5 _ K2A ]V) V-x = 

Using the following Lyapunov function, it is straightforward 
to show that this PD coupling control law drives the system 
to the desired rest state q^ globally and asymptotically while 
tending to the synchronized flow-invariant manifold: 



V 



q P 
q m 



q P 
qm 



- M1+M2 

2 

Mi-M 2 
2 

T 



Mi-M 2 
2 

Mi+M 2 
2 



q P 



(51) 



(Ki-K 2 )A 






(Kx+K 2 )A 



q P 

Vqm, 
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where Mi = M(qi), M 2 = M(q 2 ), 

q P = ^(qi + q2). 



j_ / qi + q2 
^ \qi - q 2 
The rate of V can be computed as 

^_ = _ ( qp x 7 

(it Vq m 



Ki - K 2 






Ki + K 2 



q P 



q P 
qm 



<o (52) 



which implies that V is negative semi-definite with Ki+K 2 > 
and Ki - K 2 > 0. 

By invoking LaSalle's invariant set theorem l38l . we can 
conclude that q p , q p , q m , and q m tend to zero with global 
and asymptotic convergence. This implies that qi and q 2 will 
follow q^ while qi and q 2 synchronize asymptotically from 
any initial condition. This can be extended to arbitrarily large 
networks as shown in Figure [T] 

Note that if A = 0, the PD coupling control law in (48) 
reduces to velocity coupling 



r\ = -Kiqi + K 2 q 2 , r 2 = -Kiq 2 + K 2 qi 



(53) 



This velocity coupling control can also be derived from the 
exponential tracking control law in ([5]) by setting q d = 0, 
q^ = 0, and A = 0. Therefore, the proof of the linear PD 
synchronization with A = is the same as Section IV-C 



whereas the convergence rate is now exponential compared 
with the asymptotic convergence of the PD control. On the 
other hand, we can find that positions do not synchronize in 
the absence of the gravity term even though the velocities 
synchronize exponentially fast. 

B. Synchronization with Limited Communication Bandwidth 

We now consider multiple dynamics with partially coupled 
joints (or partially coupled variables). For example, we can 
assume that only the lower joint is coupled in a two-robot 
system having two joint variables with q = (xi,x 2 ) T for 
(i = l,j = 2) or (i = 2J = l): 

Ti = M(q i )q i , r + C(qi, q*)qi,r + g(q») - Kis» 



(55) 



K 2 (^) + K 2 A ( q ) (54) 

Nevertheless, Theorems [5] and [6] are true with diagonal 
matrices, Ki, K 2 and A, which can be verified by writing 
the closed-loop system as in ([17]): 

M(q 1 )si + C(qi,qi)si + (K 1 + 
M(q 2 )s 2 + C(q 2 , q 2 )s 2 + (K x + 
u(t) = K 2 P( Sl +s 2 ), 

It is straightforward to prove that Theorems [5] and [6] still 
hold. This is because (Ki + K 2 P) and (Ki — K 2 P) are 
still uniformly positive definite, enabling exponential synchro- 
nization and exponential convergence to the desired trajectory, 
respectively. Hence, we did not break any assumption in the 
proof of Theorem [6] This partial-state coupling also works for 
the semi-contracting case with ([34]), as presented in Section 



■K 2 P)si = u(t) 
K 2 P)s 2 = u(t) 
P = diag(l,0) 



IV-D While the coupled variables synchronize to the weighted 
average of the initial conditions, the uncoupled variables 
synchronize to zero. 



C. Synchronization with Time-Delays 

The results may be extended to time-delayed couplings (see 



also |4]|). Consider for instance the two-robot dynamics (17), 
which becomes 



M(qi)si + C(qi,qi)si + (Ki 
M(q 2 )s 2 + C(q 2 ,q 2 )s 2 + (K 1 
T2i = K 2 (s 2 (t-T)-si(*)) 
7i 2 = K 2 (si(*-T)-s 2 (*)) 



K 2 )si = r 2 i 

K 2 )S 2 = Ti2 



(56) 



where the constant T > denotes the communication delay 
between the two robots, and s(t — T) = q(t - T) + Aq(t — 
T)-(q d (*) + Aq d (*)). 

Theorem 9: The robots in ( [56] ) synchronize globally asymp- 
totically VT > under the assumptions of Theorem [3] 

Proof: Similar to [40], l46ll consider the differential 
length 

V = 5z T (V T [M]V) Sz + f 5z T {e) K ° 

Jt-T 



8z(e)de 
(57) 



Using p6f shows that V is negative semi-definite, 

I/ = -2^[ k ^ Ki -kJ* z 

- 5t2 1 K.2 1 St 2 i - St^K^St^ 



(58) 



Note that Ki + K 2 > Ki - K 2 > implies K 2 > 0. 

Also, V is bounded. By Barbalat's lemma, V tends to 
zero globally asymptotically, which in turn implies Sz — » 
0, St2i — > 0, 5t\2 —> 0- Hence, the solutions converge to a 
single trajectory asymptotically. Since si(t — T) = s 2 (£ — T) 
is a particular solution of (56), qi(t) and q 2 (t) globally 
asymptotically synchronize regardless of T. ■ 

D. A Perspective on Model Reduction by Synchronization 

Another benefit of exponential synchronization is its impli- 
cation for model reduction (8). Exponential synchronization of 
multiple nonlinear dynamics allows us to reduce the dimen- 
sionality of the stability analysis of a large network. As noted 
earlier, the synchronization rate (D 2 ) is faster than the tracking 
rate (Di). Assuming that the dynamics are synchronized, the 
augmented dynamics in ( [T4] ) reduces to 

M(q)s + C(q, q)s + (K x - 2K 2 )s = (59) 

where q = qi = • • • = q p , and Di = Ki — 2K 2 is replaced 
by Ki - K 2 for p = 2. 

This implies that once components of a network are shown 
to synchronize, we can regard them as a single dynamics of 
reduced dimension, which simplifies any additional stability 
or perturbation analysis. 

VII. Simulation Results 

A. Tracking Synchronization of Four Robots 

Even though the local coupling structure of ([7]) and ([5]) has 
been emphasized, the difference from all-to-all coupling is not 
evident in a network comprised of less than four members 
(p < 3). To illustrate the effectiveness of the proposed scheme 
for a robot network with p > 4, a network of four identical 
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(a) Simulation result. 

Tracking errors of four-robot network 
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■ 1 Robotl 
Robot2 

- Robot3 
Robot4 




(b) Tracking error of each robot 
Fig. 4. Synchronization of a four-robot network 



3-DOF robots is considered here (see Fig[TJa)). The dynamics 
modeling of the 3-DOF robot is based upon the double inverted 
pendulum robot on a cart (see Fig. |4ja)), and each joint is 
assumed to be frictionless. The physical parameters of each 
robot are given in Q. 

The simulation result is presented in Figure [4] The four 
identical robots, initially at some arbitrary initial condi- 
tions, are driven to synchronize as well as to track the 
time- varying desired trajectory: Sd(t) = 0.2t, 0id(t) — 
cos(0.02tt£), 2d (t) = 7r/4(l - cos(0.08tt£)). For the control 
gains in the control law ([5]), we used Ki = 51, K2 = 21, and 
A = 51. According to Theorems [5] and [6] the tracking gain 
Ki — 2K 2 is smaller than the corresponding synchronization 
gains K1 + 2K2 and K4. Figure |4jb) shows the tracking errors 
of the four robots. We can see that the robots synchronize 
exponentially fast from arbitrary initial conditions, and this 
synchronization occurs faster than the exponential convergence 
of tracking errors. Such a result can be useful to rapidly 
achieve a collective motion of multiple robots in the presence 
of external disturbances. If Ki = 2K2 instead, we can easily 
show that the robots synchronize with global exponential 
convergence, although the tracking errors will not tend to zero. 



B. Simulation of Concurrent Synchronization for Ten Robots 

Let us consider three heterogeneous networks connected 
by feedback hierarchies shown in Fig. [5ja). If we view the 
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(a) Three networks connected by feedback hierarchies. 
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(b) Tracking error of each robot 
Fig. 5. Synchronization of ten robots on three heterogeneous networks 



reference trajectory inputs (dashed lines in the figure) as 
edges of the graph, as discussed in Section [V| this graph is 
unbalanced. The first network (robot 1 through 4) is identical 
to the network of four 3-DOF robots in the previous section. 
All the physical parameters of the second network (robot 5 
through 7) are twice larger than those of the first network 
while those of the third network are 1.5 times larger. Also, 
note that the third network is on an inline configuration by 
connecting the second feedback gain K 2 of the control law in 
([8} back to itself. All the previous theorems still hold for such 
an inline configuration (see Corollary [2]). As seen in Fig. [5jb), 
the first, second, and third network individually synchronize 
robots within each network from arbitrary initial conditions. 
The first and second network also follow the reference trajec- 
tory command from the adjacent members. For the reference 
trajectory commands, we assume that we can send the position 
and velocity values (q^, q^). Then, q^ is estimated by a high 
pass filter. Eventually, all the robots exponentially synchronize 
and then follow the desired trajectory staying together. This 
concurrent synchronization can be used to construct a complex 
and large-scale dynamic network consisting of an arbitrary 
number of heterogeneous robots and networks. 

C. Simulation of Adaptive Synchronization 
We assess the effectiveness of the proposed adaptive syn- 



chronization control law presented in Section IV-E Consider 
a two-robot network, comprised of the two-link manipulator 
robots given in page 396 of (38). The desired trajectory for 
the first joint is 0id(t) — sin7r£ and for the second joint is 
@2d(t) = 2(1 — cos0.67r£). The initial parameter estimates 
for both robots are defined as a(0) = (3, 1, 1, 1) T , and 
T = diag(0.03, 0.05, 0.1, 0.3). Figure^a) shows the synchro- 
nization of the two robots with stable tracking by Ki = 201, 
K2 = 151, and A = 101. Hence, the synchronization gain 
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Synchronization with Asymptotic Tracking 



Ki + K2 is larger than the tracking gain Ki — K2. Note 
that the synchronization of the tracking errors implies the 
synchronization of the state variables (i.e., joint 1 and joint 
2). As discussed before, the synchronization occurs faster than 
the tracking. This simulation result indicates that the proposed 
adaptive control law can be used to synchronize motions of 
robots with unknown physical parameters. 

In contrast, Figure [6jb) shows the synchronization with 
indifferent tracking by the gains of Ki = K2 = 201 and 
A = 101. The tracking dynamics then have zero gain (indif- 
ferent). So the two robots synchronize asymptotically, while 
the tracking errors remain within a finite ball. In both cases 
(Figs. [6ja) and (b)), the proposed adaptive control law ensures 
neither the asymptotic convergence nor the synchronization of 
the physical parameter estimates a^, unless the persistency of 
excitation condition is met. Nevertheless, the robots synchro- 
nize asymptotically from any initial condition in the presence 
of such parametric uncertainty. This adaptive version makes 
the proposed synchronization framework more practical. It 
should also be obvious to readers that the proposed control law 
can be straightforwardly extended to robust adaptive control 
schemes based on sliding control (38). 

VIII. Conclusions 

We have presented the new synchronization tracking control 
law that can be directly applied to cooperative control of 
multi-robot systems and oscillation synchronization in robotic 
manipulation and locomotion. We have also shown that com- 
plex dynamic networks can be constructed by exploiting 
concurrent synchronization such that multiple groups of fully 
synchronized elements coexist. The proposed decentralized 
control law, which requires only local coupling feedback for 
global exponential convergence, eliminates both the all-to-all 
coupling and the feedback of the acceleration terms, thereby 
reducing communication burdens and complexity. Further- 
more, in contrast with prior work which used simple single 
or double integrator models, the proposed method permits 
highly nonlinear systems. Providing exact nonlinear stability 
results constitutes one of the main contributions of this paper; 
global and exponential stability of the closed-loop system has 
been derived by using contraction theory. Contraction analysis, 
overcoming a local result of Lyapunov's indirect method, 
yields global results based on differential stability analysis. 
While we have focused on the mutual synchronization prob- 
lem where synchronization and trajectory tracking take place 
simultaneously, the proposed method has also been shown to 
be a generalization of the average consensus problem that does 
not address trajectory tracking. It has been emphasized that 
there exist multiple time scales in the closed-loop systems: 
the faster convergence rates represents the transient boundary 
layer dynamics of synchronization while the slower rate deter- 
mines how fast the synchronized systems track the common 
reference trajectory. Exponential synchronization with a faster 
convergence rate enables reduction of multiple dynamics into 
a simpler form, thereby simplifying the stability analysis. 

Simulation results show the effectiveness of the proposed 
control strategy. The proposed bi-directional coupling has 
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Fig. 6. Adaptive synchronization of two two -link manipulator robots 



also been generalized to permit partial- state coupling and 
uni-directional coupling. Further extensions to proportional- 
derivative coupling, adaptive synchronization, time-varying 
network topology, and concurrent synchronization of multiple 
heterogeneous networks on unbalanced graphs exemplify the 
benefit of the proposed approach based on contraction theory. 
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